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ABSTRACT 

The authors present a state feedback control approach to the single infinite bus machine and multimachine 
systems all in-cooperating a static synchronous compensator (STATCOM). The proposed controllers' designs 
are based on a linear time invariant model of the plant and state feedback scheme. First, the linear 
mathematical model of STATCOM is derived. Then, using polynomial algorithm and pole placement algorithm, 
a state feedback control law is derived. The proposed control strategy is tested on two different case studies that 
is SMIB system and three area five machine system by digital computer simulations using matlab/simulink 
program for various types of loads and/or disturbances. Comparison of these results with those methods and 
without controller establishes the elegance of these control approaches. 
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I. INTRODUCTION 

Power system oscillations happen as a result of lack of adequate damping torque at generators' rotors 
.This situation might happen as a result of a heavy loads in thepower lines, weaker interconnections, very high 
gain excitation system, sudden and drastic chane in loadse.t.c [1]. The oscillations of the rotor causes some 
oscillations from other power system variables like transmission line's active and reactive power, voltages at 
buses, bus frequency and the frequency of the oscillation is usually between 0.1 and 2Hz termed as low 
frequency [2]. There are types of oscillation: Local mode, Inter area mode, Control mode and Torsional mode of 
oscillation [3]. devices used as damping controllers for power system oscillation includes static synchronous 
compensator ( Statcom), power system stabilizers (PSS), High voltage Direct Current(HVDC) links, Thyristor 
controller series capacitor, (TCSC) etc. PSS is usually used on selected generators for damping local mode 
oscillation and sometimes used in for inter-area mode oscillation but supplementary controller is much better in 
inter-area mode oscillation damping when applied to/with the FACTS [4]. Normally the design procedures for 
these controllers are based on alinearise model under range of different operating points [4-5]. 
STATCOM is a voltage source converter based FACTS device and shunt connected similar to SVC, that are 
usually meant for voltage stability and regulation. It is also used for improving power system stability by 
exchanging reactive power to power networks [6-9]. Study of STATCOM device on voltage and power stability 
enhancement is done [10]. Liner quadratic regulator (LQR) gives optimalcontrol fora linear system based on a 
quadratic performance index minimization [11]. Many attentions havebeen put in for investigation and the use 
of state feedback based controllers by different authors such as [12-14]. 

In this research study, a supplementary state feedback based controller has been tested and proposed for 

STATCOM. Several Feedback algorithms have been utilized for a few different appropriate operating 

conditions. Two different case studies have been taken into considerations 

[1] Single machine infinite bus and 

[2] Multimachine (three area five machine system). 

The remain sections of the paper are as follows; 
Section ii discussed the statcom compensation model, relevant equations of modeling the Statcom, linescurrents 
and reactances are derived. 
Section iiioverviews the state feedback concepts 

Section iv detailed the actual power system model with a Statcom, results and discussions. Then the paper's 
conclusion was done in Section v 
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II. THEMODELING OF POWER SYSTEM WITH A STATCOM 

As seen in Fig. 1 .0, a STATCOM is connected to a transmission line via a transformer and it consists 
of a three phase gate turn-off (GTO) base voltage source converter (VSC) and a Direct Current capacitor. It's a 
reactive current source andassumed positive. 



Tl . X ll X 12 . X T2 «£ 



I 4 h cddl 



Fig.1 SMIB with STATCOM 
device 



Where 



Fig.2 Equivalent SMIB with 
STATCOM device 



/ s = (K (Vref -V )- Is) I T 



1=1+1 = (I cose + jl sin 6) (1) 

s sd sq \ s •* s / v/ 



K is again factor. 

The voltage difference between the STATCOM bus voltage, v L (t) and v 0( t) generatesactive and reactive 
power exchanges between the STATCOM and the power system, that is controlled by adjustment of the 
voltage magnitude V 0 and the phase angle' \|/' but the two voltages are in same phase so v|/ is assumed zero ( 
0°). STATCOM is installed to maintain AC bus voltage v L(t) in the power system and enhances oscillation 
damping. STATCOM control is implemented through the pulse width modulation(PWM) ratio m and phase 
angle \j/ as seen in Fig. 3 below [11], 




Fig. 3 STATCOM control technique 



dV nr I nr mk 

— = — ^- = / (2) 

^ C DC ^ DC 

From Fig.2, 

I * = I l - I .V t = JtJi + JtJz+V ( 3 ) 
Substituting eqn. I 2 into eqn. 3 gives 

f = jxj { + jx 2 (7 i -T) + v=j(x l + x 2 )7 l - jx 2 T + v (4) 

That is 

j(x 1 + x 2 )(I d +jl q ) = Vt + jx 2 T -V 

=> j(x l + x 2 )(i d + ji q ) = x i +(E' -x' i i i )-jx 2 i t (cos0 + jsin6)-v(tos5 + js\n5) After expansion and comparism of 
the real and imaginary parts, the expression of Iq and Id become;. 
E ' +x A cos 6 -V sin 8 

* d = ; (5) 

Xj + x 2 +x d 
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VcosS +x A sinO 
I = (6) 

q v 7 

X , + X n +x 

1 2 q 

v m = jx 2 i 2 +v = v + jx 2 (i d + ji q - i s ) = v md + jv mq Therefore 

y mA + i y mq = v (™ sS + i sin s ) + i x 2 K + i 7 ? - 7 S ( C0S * + i sin *)} ( 7 ) 
. Substitution ofequation. 5 and equation 6 into equationalso comparing the real and imaginary parts gives us 

(x + x' )V cos 5+ E ' x , + / cos^jc^fjc+jc' ) 

y mq = — — — *- J — L — — — (8) 

x 1+ x' d +x 2 
(x,+ x )V sinS+ I sinO x Jx ,+x ) 

y md = 2 -^- (9) 

X , +X + X ^ 

1 q 2 

E'V v 2 x'-x 

P e = q m sin 0 + — sin 26 (10) 

Xj+ x' d 2 (Xj+ x q )(x 1 + x' d ) 

The dynamics of the generator and the excitation system are expressed as a fourth order model and given as 

AS = co b Aco (11) 

Acq = -(AP + DA cd) I M (12) 
*K = (-AE' q + (x d -x' d )Ai d + AE fd )ir; Q (13) 

AE fd =-^AV t --j-AE fd (14) 

/ s = (K Au - AIs) IT (15) 

Let the input of STATCOM controller to be;. 

Au = {V ref -K u AV m +K w Aco) (16) 

^ M and ^are the gains of voltage and damping control loop, respectively, V re fis the reference voltage of the 

STATCOM regulator andKr Gain of the stabilizing signal. 

Thus 

Au = -K r K 4 K u AS + K r K a Aa - K r K 5 K u AE' q - KK 6 K u AI s (17) 

Where K } -K i5 are linearization constant and brings the system toalinearised power system model in matrix 
form (18) below 
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( (18) 



III. STATE FEEDBACK CONTROLLER 

Pole placement with a state feedback 

Assume that the single-input system dynamics are given by 
Table 1.0 Butter worth Polynomial factors 



(19) 



x(t) = Ax(t) + Bu(t), 

y(t) = Cx(t)+ Du(t). 
Assume a control law of the form: 

u(t)=-Fx(t) or u(t)=r(t)-Fx(t) (20) 
Is called state feedback and brings system to the form [15] 
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x(t) = (A - BF)x(t) + Br(t), 

(21) 

y(t) = (C - DF)x(t) + Dr(t) v ; 

Closed loop transfer function 



0— (*4 ,J 



3 



Fig 4.The closed loop control system 

We have: 

x = (SI-A y 1 B(r-Fx) o sx-Ax= Br-BFx (22) 

Therefore, y= Cx, now closed loop transfer function becomes;. 

s d (s) =C (sI-A F ) *B where A p = A-BF ,(23) 
And the Closed loop characteristic polynomial is: 

x c =det(sI-A F ) (24) 

In brief 

Assume(A,B) are controllable, then for any arbitrary polynomial 



x ,(s) = s + a , s 

cl v ' n-1 



.+ a , s + a , 



(25) 



There willexist a state feedback gain K, that x cl (s) is closed loop characteristic polynomial, i.e., 
x cl =det(sI-A p ) (26) 

The matrix A F = (A-BK) reveals not only the stability of the closed loop system, but also a response of the 
rotor angle or speed output. Thus, the controller gain, K, matrixmust be designed a such that will meet the 
stability criteria which will be achieved by means of placing the closed loop poles at a desired locations in 
the left half of a complex plane. K can therefore be determined using ng: 

(1) Equating of equation coefficients 

(2) Pole placement algorithms 

In this paper the input to STATCOM's controller is taken as Vref and the outputs are taken as from the output 
of the modeledpower system (rotor speeds, real and reactive powers, rotor angles etc.) eachunder various 
operating conditions. 

2. Polynomial Placement Algorithm 

It is from an evidence [15] that with somehigher order systems, ordinary pole placement can become tedious in 
dealing with, one method for placing poles is by employing the polynomial approximation technique such 
asITAEand Butterworth pattern. 

In these ApproximationsAlgorithms procedure is summarized as: - 

• Determining the desired settling time ts - 

• Find k = n polynomial from the table 1 or 2. - 

• Divide the pole locations by ts - 

• Form desired characteristic polynomial Od(s) 

• andapply acker/placematlab commands to know the feedback gains. - 
Simulate to study the performance and control efforts. 
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Table 1.0 BUTTERWORTH Polynomial factors 



n 


Denominator coefficients for polynomials 


1 


(s + 1) 


2 


( s 2 + 1.4142s + 1) 


3 


(s + l)(s 2 + s + 1) 


4 


(s 2 + 0.7654s + l)(s 2 + 1.8478s + 1) 


c 

J 


(S+1 ){S + U.OloUS + 1 ){S + l.OloUS + 1 ) 


6 


(s 2 + 0.5176s + l)(s 2 + 1.4142s + l)(s 2 + 1.9319) 


7 


(s + l)(s 2 + 0.4450s + l)(s 2 + 1.2470s + l)(s 2 + 
1.8019s + 1) 



Table 2.0 ITAEPolynomial factors 

1 5 

s + co 0 

S 2 + 1 .4 03 Q S + 03 2 

s 3 + 1.7 5o3 0 s 2 + 2.15o3 2 s + 03 3 

s 4 + 2.1o3 0 s 3 + 3 .4 03 2 s 2 + 2 .7 03 3 s + 03 4 

o5 „ „ 4 .-2 3 ,_ ,_ 3 2 ^ . 4 1 5 

S +2.8o3 0 s +5o3 0 s +5.5o3 0 s + 3 .4 03 0 s + 03 Q 
S 6 + 3.2 5 03 0 s 5 + 6.6 03 2 s 4 + 8.6 03 3 s 3 + 7.4 5 03 4 s 2 
+ 3.9 5o3qS + 03q 

S 7 + 4.4 75 03 0 s 6 + 1 0.4 2 03 2 s 5 + 1 5.0 8 03 3 s 4 + 1 5.5 4 03 0 V 
+ 10.6 4o3qS 2 + 4.5 8o3qS + 03^ 

A. Pole placement algorithm 

As K=fV T the elements of F are selected while the process of finding V is as follows [16]: 

n n-\ n-2 

Let the characteristics equation of the system matrix be* + a A s + a 2 s + a n 

For closed loop purpose ,it is required to move the eigen value of the uncompensated system characteristics 
equation as soon as eigen value of the closed system 7 1 > 7 2 , 7 1 > 7 2 , V 3 V n are specified the desired 

closed loop characteristics equation is found from 

// 

PU)=M {1-7, )=(s-7,){s-7 2 ) (s-7 a ) (28) 

i=l 

The difference between the open-loop Xtics polynomial and that of the desired closed loop Xtics polynomial is 
D. 

D(s)=P(s)-A(s) (29) 
D(s)={p-a 1 )s n -' +(p-a 2 )s n - 2 + (p-a n ) (30) 

Then we defined a vector d whose element the coefficient are of d ( s ) 

* = [(p, -Mp, ■«.) {p.".)] (3D 

Select a vector F such that the element of the matrix 

Q = [bF;ABF;A 2 BF A" 1 B F ] (32) 

Is completely state controllable, then the vector V can be found as 
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-1 — 1 T 

V = Q X d 



Where X is a Toeplitz matrix defined as 



(33) 



r i 



X = 



o 
1 



o 0~\ 

o o 
o o 



'J 



(34) 



B. PID-STATCOM Pole Placement Approach 

In the design of the PID controller, the gain settings can be computed by placing the Eigen values at a 
pre-specified locations, this is usually known as the pole placement method where K p , Ki, K D , are the gains of 
the PID controller and T w is the wash out time constant as shown in fig 3, 

The approach starts with linearizing the non -linear model around a nominal point to obtain the desired linearized 
model described by input-output equations 












— ► 


s 







Fig. 5 The magnitude control block diagram 




Fig. 6 PID-STATCOM Based power system closed loop block diagram 
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= Ax(0+Bu{t), (35) 
y&) = Cx(t) + (36) 

Where x, y and u are the state vector variables, output signal and control signal respectively, taking the Laplace 
transform and substituting the output equations into the state equation gives 

X(s} = (sf - A^BlKs) (37) 

If the system output Y is fed into the input of the PID controller as shown in fig 4, and the PID controller is 
having a transfer function H(s), we can therefore write the control signal as 

U(s)=H(s)Y(s) where 

H ^ = ^ ^P+Y+W (38) 
Putting (38) into (37), we have: 

Xis} = (si - Ay L BH(s} CX(s}(39) 

Or 

[1 - (si - Ay L BHls}C] X(s} =0 (40) 
If X is the assigned Eigen value of the whole closed loop system, then 

Det[l - QJ - Ay L BH0^C=0 (41) 
Equ.15 can be re- written as 

l-CQJ -Ay L BH03=0 (42) 

or 

A (*) = — 7Z ( 43 ) 

Equ. 38=Equ.43, thus 

" w k + Tw k + -^^-k.q = r44^ 




4.1 Three Area Five machine system, with STATCOM 
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IV. RESULTS AND DISCUSSION 

The effectiveness of the proposed methods was tested on two case studies 

(1) Single machine infinite bus system. 

(2) Three Area five machine system 
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Reactive Power for line 21-25 
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Fig. 8 Reactive Power for m=line 1 lto 5 



,-3 Speed deviation of various controller for case 2 
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Fig. 9 Speed deviation of case2 



Angle response for various controller easel 
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Fig. 10 angle respone of case 12 
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To assess the effectiveness of the proposed methods, two different loading conditions are considered 
which name easel and case 2. The results for the systems are presented in what follows: 
Case 1 (normal loading): Operating points of SIMB= Pe=1.0pu at unity pf.; Vt=1.0 p.u. 
Case 2 (heavy loading): Operating points of SIMB= Pe=1.2pu at 0.8 pf.; Vt=1.2 pu. 

multi-machine power systems stability enhancement by means of ST ATCOM -based stabilizers was studied. To 
demonstrate the influence of the introduction of STATCOM [28] to the multi machine power systems,The three- 
area five machine system model shown in Fig. 8, is used, the STATCOM is inserted between buses 21 and 25. 
The voltage at bus 21 and reactive power flow into bus 25 is controlled is controlled by the statcom element. 
Generators 1 to 4 have static exciters with power system stabilizers, and generator 5 has a dc exciter. Each 
generator has a thermal turbine with governor. There is one under damped low frequency electro -mechanical 
mode in which generator 5 oscillates against generators 1 to 4 

A comparison of the responses obtained by the three proposed method based on pole placement techniques, 

PID controller tuning based on pole-placement technique produce optimum controller functions for linear 

systems designed for case 1 and case 2 operating points. By assigning three poles of the compensated closed 

loop system the appropriate values of K P , Ki and K D has been obtained. The location of the dominant 

eigenvalues were selected to be -1.8561 ± j8.2953, corresponding to the damping ratio of 0.2186, the PID 

controller parameter for the nominal operating point(case 1) is obtained to be: KP=-4.95; KI=-296.3 and 

KD= 113.23. While for pole algorithm described in section 3 the controller parameter are obtained as [Kl, K2, 

K3, K4, K5] = [-0.9680 63.7854 -27.6727 -0.1683 0.1585]. For ITAE 

While for ts=5s , we have -1.6054, -0.8672 ±j 1.7504, -1.3143±/ 1.1357, -1.5365±/ 0.5616i 

The ITAE controller parameter are obtained as [K b K 2 , K 3 , K 4 , K 5 , K6, K7] = [32.9, -9.2,420.9,-18.2,2660,- 

4058,-1066]. Figure 10 shows the response of real power in line 21-25 where statcom is injecting reactive power 

from bus 25 for different algorithm based on pole placement. Figure 12-13 shows the response for reactive for 

ITAE, Butterworth and PID 

controller. Figure 1 1 demonstrates the voltage response of bus 25 where the statcom is connected while figure 
12 shows the speed deviation of (G1-G3) response under different control algorithm. 

V. CONCLUSIONS 

In this paper, STATCOM Controller based on state feedback concepts are proposed for damping 
oscillations and the effectiveness of the proposed control methods are compared within themselves under some 
disturbances. The controllers are tested on single machine infinite bus System From the The 

desired poles for 7 th order reduction are-8.0271, - 4.3361±/8.7519, - 6.5714±/5.6786, - 7.6824 ±/2.8081 

results it can be concluded that the state feedback based on ITAE produces no steady state error and acceptable 
overshoot under some disturbances. 
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